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ABSTRACT 
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This work examines the discrete spectrum of the Orr-Sommer feld problem 
of hydrodynamic stability for boundary layer flows in serai-infinite regions. 
Related questions concerning the continuous spectrum are also addressed. 
Emphasis is placed on the stability problem for the Blasius boundary layer 
profile, 

A general theoretical result is given which proves that the discrete 
spectrum of the Orr-Sommerfeld problem for boundary layer profiles (U(y), 
0,0) has only a finite nimber of discrete modes when U(y) has derivatives 
of all orders. This has been suspected for some time on the basis of numer- 
ical evidence, but the present result is the first theoretical proof of its 
type for unbounded flows. 

Details are given of a highly accurate nimerical technique based on 
collocation with splines for the calculation of stability characteristics. 
The technique includes replacement of "outer" boundary conditions by asymp- 
totic forms based on the proper large parameter in the stability problem. 
Implementation of the asymptotic boundary conditions is such that there is 
no need to make apriori distinctions between subcases of the discrete spec- 
trum or between the discrete and continuous spectrums. Typical calculati^fus 
fo.: the usual Blasius problem are presented. Results also show that there 
are not a large number of discrete temporal modes of the Blasius problem 
lying close to the temporal continuun. 

The parallel flow assumption, which leads to the usual Orr*Sommerfeld 
problem, ignores the small V velocity component of the mean profile. If 
V is retained, a modified Orr-Somraerfeld equation involving first and third 
derivative terms is obtained. This modified equation has recently proved 
competitive with non-linear approaches in predicting instability for a model 
profile where V is stabilizing. Ihe modified Orr-Sommerfeld problem is 
considered here for the Blasius problem in which V is destabilizing. Both 
the marginal stability curve and the higher modes of the modified problem 
are discussed. Critical para&eters for instability are lowered by 3Z, but 
this is not enough to reconcile the linear theory with experimental results. 
A conjecture is given as to when linear stability predictions with V 
included do not need further refinement from multiple scale methods. 
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INTRODUCTION 
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In hydrodynamic stability, a key assumption in the linear theory for 
the stability of small amplitude disturbances is that the disturbance stream 
function can be expanded in terras of a complete set of normal modes. These 
modes, in turn, are obtained from the eigenfunctions and eigenvalues of an 
appropriate boundary value problem. For example, for boundary layer type 
flows over a flat plate, the parallel flow assumption and consideration of 
two-dimensional disturbances leads to the usual Orr-Soimner feld (0-S) problem 
for To llmien-Schl itching waves. This problem involves a fourth-order linear 
differential equation and boundary conditions which come from the no-slip 
requirement on solid walls and appropriate conditions at infinity in the 
free stream. 

Several basic questions arise naturally in the study of linear 
stability problems. The first question, of course, concerns the stability 
of the mean profile being considered. This is really a question about the 
first, or least stable, eigenvalue. As, in general, no closed- form 
solutions of the eigenvalue problem are possible, analytical work exploits 
the fact that the Reynolds number R is large for most flows of physical 
interest. Hence perturbation methods can be used to obtain asymptotic 
approximations to the eigenvalue relation, and a curve of marginal stability 
can be derived. On this curve, the disturbance wave number oi, wave speed 
c, and frequency o) - ac are all real, and hence the temporal and special 
instability problems coincide. 

Inside the marginal stability curve, small amplitude disturbances will 
grow either specially (a complex, o) real) or temporally (a real, c com- 
plex). Numerical methods can be used to obtain curves of constant amplifi- 
cation and constant phase speed in the (a,R) plane for the temporal problem 
or the (m,R) plane for the special problem. Most work has focused on the 
temporal case where the eigenvalue c enters in a simple way. However, for 
comparison with experimental measurements, the relatively complicated spe- 
cial case is probably more appropriate. 

Results for marginal stability and curves of constant amplification and 
phase speed are of great importance, both theoretically and for application 
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to the practical engineering problem of transition prediction (Obremski, 
Morkovin, and Landahl, 1969). However, information of this sort is not 
adequate to attack the problem of boundary layer »-eceptivity . Consider, for 
example, the Orr-Sonmerfeld problem for the usual Blasius boundary layer 
profile on a flat plate ignoring the outflow velocity from the bound- 

ary layer. The receptivity problem requires the proper decomposition of an 
external disturbance into a superposition of linear stability modes. 

In particular, to use linear stability theory to study the evolution of a 
disturbance, we must know how the initial amplitudes of modal disturbances 
in the flow are caused by external disturbances imposed on the flow. 

Although few studies consider this basic question, the need to address it 
has been realized for some time. Mack (1977), in a discussion of the appli- 
cation of stability theory to transition prediction for boundary layers 
states that **. . . if there were no disturbances, there would be no transi- 
tion and the boundary layer would remain laninar. Consequently, it is fu- 
tile to talk about transition without in some way bringing in the distur- 
bances which cause it . . . .** Further on. Mack states that **. . . the pre- 
cise mechanism by which, say, free stream turbulence, sound and different 
types of roughness cause transition remains to be discovered." 

Resolution of this type of question involves a form of wave packet 
analysis. The general problem is the solution of the initial- boundary value 
stability problmn in which a disturbance ^ich is initially localized in 
space and time interacts with the mean flow, propagates in space, and 
evolves in time. The theoretical question is to determine the distribution 
of energy from an arbitrary disturbance among the of the Orr- 

Sonmerfeld equation and to calculate the subsequent evolution of the 
disturbance* 

In theory the solution to this problen is quite straightforward. Vfe 
simply determine the complete set of normal modes of the Orr-Sommerfeld 
equation and expand the initial disturbance in terms of them. The theory 
required to formally carry out this expansion has been recently derived by 
Salwen (1980). Once the disturbance has been expanded in terms of the nor- 
mal modes, the initial amplitude of each mode is known and the evolution of 
the disturbance wave packet in the flow is completely determined by the 
evolution of each mode. 



To carry out the expansion of an arbitrary disturbance in terras of 
normal modes of the Or r*-Sommer feld equation requires knowledge of the full 
spectrum for the boundary layer flows on unbounded domains. Unfortunately, 
key information on modes in the discrete spectrum has been lacking in both 
the temporal and special problems. It is, of course, well known that there 
is at least one discrete mode because a neutral stability curve exists. 
Grosch and Salwen (1978a, b) have shown that both the temporal and special 
stability problons in an unbounded domain always have a continuous spectrum. 
What is not known, and is vital in order to carry out the expansion, is how 
many discrete modes exist. This report contains a theoretical proof that 
there are only a finite number of discrete modes. Ihe available evidence, 
nimerical, further suggests that the nimber of discrete modes is relatively 
small, see, for example, Jordinson, 1970; Mack, 1976; Corner, Houston and 
Ross, 1976; Murdock and Stewartson, 1977. All of these investigators, and 
others, find a large number of possible discrete modes near the continuum. 
Most authors, however, interpret these modes as spurious discrete modes. 

The major exception is the work of Antar and Benek (1978) which claims to 
have found two families of discrete modes lying close to the temporal 
continuum . 

To distinguish between a valid discrete eigenvalue and a continuum 
eigenvalue in either the spacial or temporal case requires accurate determi- 
nation of the associated eigenfunction from the boundary out to a large 
distance into the free stream. If y denotes distance from the solid lower 
boundary, the anplitude ^(y) of the disturbance stream function for a 
discrete mode will be zero at y * 0, reach a maximum inside the boundary 
layer, and decay exponentially as y ^ By contrast, a continuun eigen- 
function is small in the boundary layer region, relatively large with oscil- 
latory behavior above the boundary layer, and merely bounded as y ^ For 
various reasons, both shooting and expansion techniques experience difficul- 
ties in computing eigenfunctions on unbounded domains. Consequently, it has 
not been clear that previously used numerical techniques could distinguish 
between discrete modes on or very near the continuum and continuum func- 
tions, so a large number of discrete modes near the continuum has remained a 
disputed possibility. The present work resolves this controversy by ruling 
out discrete modes near the continuum. The numerical technique employed 
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here is local in nature, does not produce spurious modes, and allows accu- 
rate determination of eigenfunctions in both the spacial and temporal cases. 
Results are given in Section 4. 

The Orr-Sommer feld problem occupies a central place in the linear the- 
ory of hydrod 3 mamic stability and has provided great insights into the 
nature of instabilities. However, it is important to note that the usual 
Orr-Sommer feld situation is somewhat restrictive, even within the context of 
linear stability theory. This is due to the parallel flow assumption on the 
mean velocity profile in the derivation of the usual Or r-Sommerfeld equa- 
tion. Under this assunption, the 0(R"^) mean outflow velocity is neglected 
entirely compared to the 0(1) mean streamwise flow. This amounts to 
ignoring all x-derivat ives of the mean flow streao function t(x,y,t). In an 
attempt to bridge the gap between the minimum critical laynolds minber 
for instability predicted by the Orr-Sommerfeld problan and the lower value 
of Rc found experimentally, some authors e.g. Saric and Hayfeh (1975), 
have retained full x dependence in derivatives of ^ and considered the 
stability problem using multiple scale techniques. There is, however, a 
middle course which remains within the franework of linear stability theory 
yet still partially relaxes the parallel flow assumption. Non-parallel 
effects due to outflow from the boundary layer region may be included in 
linear theory by retaining Bf/dx and neglecting only second and higher 
order x-derivatives of The accuracy of the mean flow representation in 

the stability problem now becomes the same as that of the Prandtl boundary 
layer equation. The resulting linear stability problem now involves a 
modified Orr-Sommerfeld equation with important first and third derivative 
terms attributable to the outflow velocities. For a test problem involving 
the asymptotic suction boundary layer profile, Lakin and Reid (1982) showed 
that including the 0(R~^) suction component in the mean profile signifi- 
cantly increased the minimmi critical Reynolds nunber for instability. In 
the more realistic stability problem for the Blasius boundary layer, the 
effect of the outflow component of the mean profile should be destabilizing. 
Barry and Ross (1970) have exmnined the change in Rc induced by including 
the mean outflow component. However, very little work has been done on 
higher modes of the modified Orr-Sommerfeld equation for the Blasius problem 
including outflow. Corner, Houston, and Ross (1976) consider the higher 
spacial modes, but their results for the modified Orr-Sonmerfeld problem are 
limited and inconclusive. 



The highly accurate nvjmerical procedure developed during the course of 
the present work has been used to re-examine the marginal stability curve 
for the Blasius problem including outflow. Results were also obtained for 
the higher modes of this problem in both the temporal and special cases. 
Details are given in Section 5 of this report. 

THE NATURE OF THE DISCRETE TEMPORAL SPECTRUM 


In order to use eigenfunction expansion techniques in stability 
calculations it is necessary to know the number and distribution of the 
eigenvalues of the Orr-Soranerfeld equation. In addition to the continuous 
spectrun (Grosch & Salwen, 1978; Salwen & Grosch, 1981), there is, of 
course, a discrete spectrum. The work of Jordinson (1971), Mack (1976), and 
Murdock & Stewartson (1977) suggests that the number of discrete eigenmodes 
is finite and, at finite Reynolds nuaber, small. 

Some light can be thrown on the question of the number of eigenvalues 
of the Orr-Sommerfeld equation by using a technique developed by Lidski and 
Sadovnick (1968), They considered the eigenvalue problem for the Orr- 
Sommerfeld problem in a finite domain and developed formulae for the suns of 
integer powers of the eigenvalues of the stability problem, Ihey suggest 
that these formulae could be used to calculate the eigenvalues. It seems 
that this approach to the nunerical problem is not very fruitful. However, 
their methodology can be used to address the question of the number of 
eigenvalues of the stability problan in the infinite domain. 

Here we adapt the method of Lidski & Sadovnick, (hereafter L & S), with 
some slight change in notation. The Orr-Sommerfeld problem is, for the 
discrete modes, 

- ioR{(U-C)L - U”}*, (1) 

♦ (0) • ♦ '(0) ■ 0, ♦ ♦ ♦ ' ♦ 0 as y ♦ • (2a-d) 

•►U T ^ 2 

with L “ — - 


( 3 ) 



and U(y) the velocity profile. 
Following L & S, define 


cr'" 

Or i 


I 
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r(y) = - ia R(U(y )-l) , 


(A) 


q(y) H ioRU", 


(5) 


2 ? “ iaR(c-l). 


( 6 ) 


'few L & S show that if r(y) and q(y) have derivatives of all orders, 
then the four linearly independent solutions of equation (1) are given, 
formally, by the four series: 


I r'' ^(y) 

k-0 ^ 


(7 a) 


♦2 - I (-D^-^Cy). 

k-0 


(7 b) 


-k. 


k-0 ’ 


(7 c) 


♦4 I 


k-0 


(7d) 


with z taken to be the root of (6) with positive real part. 

The functions {a (y)} are the solutions of the recurrence 

• • ^ d 

differential equations (primes denoting ly’ 
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ORIGINAL PAGE IS 

a;; » - (j) a|;_l + y (a^ - r) OF POOR QUALITY 


* ’t"' ^-3 *7 

(8) 

with a, (y) =0 if k < 0, 
k 

(9) 

ao (0) - 1; aj^(O) - 0, k > 1. 

(10a,b) 

The first few of the ^ can be found quite easily, 

k ■ 0: ad(y) ■ 0, ao (0) * 1, 

We have, for 

(lla.b) 

ao(y) “ 1; 

(12) 

k ■ 1: ai' (y) ■ y (o^ - r(y)), ai (0) ■ 0, 

(13a,b) 

, y 

«1 (y) • y ^®^y "I r(5 )dC); 

(14) 

k » 2: 


s 1 y 

(y) ■ (^)r' (y) ♦ - r(y))(o^y r(E)dC), 

(15a) 

a2 (0) ■ 0, 

(15 b) 


y 

%(y) - (^(r(y) - r(0)) + y -V ~ ^ 

V C y C 

♦ |a2 I [ I r(C)dC]« +y| r(C)[ | r(C)<lC]dC. (16) 
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The functions {b .(y)} are the solutions of the recurrence 
^ f J 

differential equations 

» [l 2 ♦ r(y) L + q(y)]b^_, 


k-2, j 


(17) 


9 

i 


with b . = 0 if k < 0, 

k,J 


and bo 1 (0) ■ 2 (0) ■ 1, 




ORIGINAL PASS ! i’ 

OF POOR Q'JALI'rY 


(18) 

(19a, b) 
(20a, b) 


Again, Che first few of the (b ,(y)} can be found. We have, for 

iC 


(2 1 ) 


k ■ 0, 


Lb . ■ 0 
0 .J 


( 22 ) 


and Che boundary conditions (19a,b). Thus, 


bo 1 “ e , 

% 


(23a) 


ho 2 ■ ® « 


Oy 


(24a) 


From (21) 


b, , 5 0. 
* .J 


(25) 


Finally, for k ■ 2 


Lb, ,(y) - [l^ ♦ rL ♦ q]bo 

J 1 J 


with b 2 .(0) - b^' .(0) - 0. 
> J tj 


(26) 
(27a, b) 
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Therefore 


OKIGINAl. FAGt \S 
OF POOR QUALFTY 


9 

y 


. (2c.)-' I [e (28 

t 

b2 2(y) - (2 o2)-' ] '->'>1 e“'S( £) SC. (29 

The structure of the (a (y)} and the {b .(y)} is now clear. Of the 

^ > J 

four solutions f'I'jj}, * I,... ,'4, we see that ((ij and .^4 grow 
exponent ionally as y + «• and ^2* ^3 decay exponentially as y > «». 

Thus the solution to the Orr-Sommerfeld eigenvalue problem is a linear 
combination of 1^2 and ^3. 

The eigenvalue relation for c is therefore 


^ 2 ^ 0 ) 

♦J( 0 ) 


•>3(0) 

♦^( 0 ) 


- 0 . 


(30) 


We have 


♦ 2 ( 0 ) ■ e *^{ao(y)- ai(y)/2 ♦ a 2 (y)/*^ ♦ 0(z~3)}, (31) 

“ **o i^y^ ^2 i(y)/*^ * o(z“**). (32) 

> » 

Using the definitions of the {aj^(y)} and {bj^ l(vy} it is easy to see that 
♦ 2 ( 0 ) = ♦ 3 ( 0 ) = I, (31) 

♦i(0) • - z ♦ {(<*2 - r(0))/2z ♦ 5 r'(0)/4z2 ♦ 0 (z'^)}, (32) 

♦ |( 0 ) S - a. 
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Therefore the eigenvalue relation (30) ia ^ TOOT> 

-a * z - {(o2 - r(0))/2z + 5 r’(0)/4z2} + 0(z“^)} - 0 

For |z| large, equation (33) is an asymptotic approximation tc the 
eigenvalue relation. From (6), |z| may be large in two ways: either aR 
large or o -nd R fixed and |c| large. Once z is obtained by solving 
(33) to the appropriate crder, c is easily obtained from (6). For 
exaaple, retaining terms in l23) to 0(z”M gives 

c ■ 1 - ia/R (34) 

which is a single point on the temporal continuua. This is not a contra* 
diction as the outer conditions (2c, d) that 0 as y > ** automatically 

imply that 4 and remain bounded as y * and these wre outer conditions 
for the continuua. To the next order, (33) is 

- az - Y (■? ~ ioR) ■ 0. (35) 

The condition Re(z) > 0 now produces a single mode «rith |c| ■■ 0(1) idiich 
does not lie on the continuum. The key point here is that equation (33) 
does not involve oscillatory or other type functions which might produce an 
infinite sequence of zeros, and hence, at each order, (33) will give only a 
finite number of modes. Further, if large |z| is identified with large 
|c I , then the fact that (33) to 0(z~^) gives (34) shows that there are no 
discrete modes for fixed a and R outside of a circle with 0(1) radius 
centered at the origin in the complex emplane. The eigenvalue relation is 
an entire function of *,e, and R and consequently, in the finite region 
of the c-plane inside the above circle, there can be only a finite nusber of 
discrete modes. Therefore, although this argusent does not give the nusber 
of discrete swdes, it does prove that the Orr-Sommerfeld problem for bound- 
ary layer profiles U(y) (which have derivativt^ of all orders) has only a 
finite ntnber of discrete temporal nodes. 



NIMERICAL METHODOLOGY AND ASYMPTOTIC BOUNDARY CONDITIONS 


3 
i 

There have previously been two basic methods for nunerically solving 
Che Orr-Sommerfeld equation; the shooting method (Mack, 1965, 1976) and the 
expansion method (Grosch and Salwen, 1968; Orszag, 1971; Salwen and Grosch, 

1972). The shooting method is quite efficient if the approximate location 
of the eigenvalue is knovn. There are, however, problems in calculating the 
eigenfunctions with this method because of the necessity of purifying or 
orthonormalizing the separate pieces of the eigenfunction. The higher 
discrete modes and the continuum eigenfunctions are especially diffic’^lt to 
calculate because they are rapidly oscillating functions of the distance 
frcm the boundary outside of the boundary layer and rapidly decaying inside 
the boundary layer. Shooting methods also explicitly require the imposition 
of the particular outer boundary conditions for either the discrete modes or 
the continuum functions. Finally the integrations must be iterated until 
the correct eigenvalue (the frequency in the temporal stability problem or 
the wave-nunber in the spatial stability problem) is obtained. 

The expansion method presents problems of a different kind. The 
solution to the Orr-Sommerfeld equ,.tion is expanded in some set of functions 
and the expansion is truncated after a finite nmber of terms. This reduces 
the problem to finding the eigenvalues of a matrix. Once the eigenvalues 
are determined it is, at least in principle, straightforward to compute the 
matrix eigenvectors and the corresponding eigenfunctions of the Orr- 
Sommerfeld equations. 

This method requires that the expansion functions be complete; and this 
requires, of course, the specification of some boundary conditions at the 
"outer" boundary. There are two ways to handle a boundary at in .nity: 
either to apply the boundary conditions at some finite distance from the 
boundary where the base flow is approximately that at infinity or to map 
the Infinite region into a finite region (Grosch and Orszag, 1977). In 
either case the boundary conditions for only the discrete modes or for only 
the continuum functions must be applied. 
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Then N eigenvalues of the N by N matrix are readily found, say by 
using the QR algorithm. Some of these N eigenvalues are approximations 
to the true discrete eigenvalues of the Orr-Somraerfeld equation and, if 
there are only a finite number of discrete modes, M, and M < N, the 
remaining modes are spurious. 

The numerical procedure used in the present wrk utilizes a solution 
process involving collocation with B-splines. The basis of this technique 
is the collocation code COLSYS developed by Ascher, Christiansen, and 
Russell (1978) for linear and non-linear boundary value problems. Because 
of the nature of the solution process, COLSYS does not suffer from the 
purification, normalization, and spurious mode problems which have plagued 
current shooting and expansion techniques. As implemented in the present 
study, there is also no need to treat discrete and continuum cases 
differently ii. terms of conditions at the outer boundary. 


As developed by Ascher et al., COLSYS is a real variable code on a 
finite interval for problems not involving eigenvalue determination. By 
contrast, the 0-S problem for boundary layer profiles is an eigenvalue 
problem involving complex variables on a semi- infinite interval. These 
differences help to explain why this collocation code has not previously 
been used in the stability context. To overcome these apparent limitations 
and produce a package which gives highly accurate nunerical results for both 
eigenfunctions and eigenvalues in the boundary layer stability context, 
several strategies were used. 


To enable the basic code to be used for eigenvalue problems, the 
eigenvalue was considered to be a varible and a scalar, first order equation 
da 


.dc ^ 

W 


or ^ * 0) was appended to the usual Orr-Sommerfeld problem. 


This changes the formulation from a linear to a non-linear boundary value 
problem, but introduces no additional difficulties so far as COLSYS is 
concerned. Indeed, increasing the order of the problem allows specification 
of an additional condition which provides a natural vehicle for normalizing 
the nunerical solutions. 
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The package of subroutines which comprise COLSYS involves only real 
arithmetic, and, because of the nature of these routines, COLSYS itself 
cannot be re-written to employ direct complex arithmetic. The complex- 
valued equations of the stability problem could be rewritten in real-valued 
form by explicitly deriving real and imaginary parts, but this procedure was 
rejected as being cunbersome, inflexible, and subject to manipulative 
errors. Instead, COLSYS itself was left untouched but the manner in which 
complex numbers are stored in computer memory was exploited so that complex 
arithmetic could be used directly in user-supplied subroutines specifying 
the eigenvalue problem. Real and complex derivatives were related using 
complex-valued Jacobians. Compatibility of COLSYS with the complex-arithme- 
tic subroutines also required appropriate re-ordering of the dependent vari- 
ables and their derivatives. 

The principal difficulty in adapting COLSYS to boundary layer stability 
involves reduction of the computational domain from the semi- infinite 
physical interval [ 0, to the interval [o,y] where Y is not excess ivly 

large, e.g. Y < 20 say. Simply putting a ”top** on the problem is known to 
produce spurious modes in the discrete spectrum. The semi-infinite interval 
may also be mapped into a finite interval [Crosch and Orszag (1977)] but the 
resulting problem is quite complicated. The present work employs asymptotic 
outer boundc.ry conditions which are equivalent to an asynptotic-mmerical 
matching at Y of the computed eigenfunction for y < Y a 'd the proper 
asynptotic form of the eigenfunction for y > Y. 

The heart of the asymptotic boundary conditions is the theoretical 
result of Lakin and Reid (1982) that unbounded domain effects induce higher 
order corrections to solutions of the eigenvalue relation for bounded values 
of y, right down to y ^ 0. Inclusion of these effects requires 
derivation and use of a proper large paraneter which differs from the usual 
large parameter in the bounded domain case [Lakin and Reid (1970), Lakin, 

Ng, and Reid (1978)]. To obtain the proper parameter, consider the Orr- 
Sommerfeld equation in the explicit form 


(E? - o2)2^ . ioR{ (0-C)(D^ - o2)<^ - U"*} 


(3.1) 



where again U = U(y) is the mean velocity profile, « is the wavenumber, 
c is the wave speed, <()(y) is the amplitude of the disturbance wave 

function, and D ~ 'fhe Reynolds number R in (3,1) is based on the 

free stream velocity Uq and the length scale L rather than the boundary 
layer thickness d ,* Length scale conversions can be made, if desired, 
using the relation 

6* = 1.7207L. (3.2) 

In the present scaling, as y U(y) 1, U**(y) 0, and the 0-S 

equation tends to a constant coefficient equation which can be written in 
the factored form 

(D^ -a2)(i)2 -X2)^ = 0 (3.3) 


tAiere 


» iaR(l-c) + and Re(X ) > 0. (3.4) 

Solutions of equation (3.3) which do not violate the outer boundary 
conditions are e"'*^ and e For continuun modes, ^ich have Re(A) = 

0, the solution e^^^ is also admissible. The e solution can be 
identified with analytic continuations of solutions of the Rayleigh equation 

(U-c) (d2 - 0 (3.5) 

to large y. As (3.5) is obtained from (3.1) by formally letting R ♦ 

^Xy 

these solutions exhibit inviscid-type behavior. By contrast, the e solu- 
tions have viscous-type behavior. Tliis argument also indicates that the 
proper paraoeter in the usual 0-S problem for boundary layer flows is A^ 
rather than i * R. In particular, the scaling in A is relative to the free 
stream velocity rather than the shear in the critical layer as in the case 
of bounded channel flows. 

Asymptotic outer boundary conditions are based on solution' of equation 
(3.3) and avoid the spurious mode difficulties known to occur in the case of 



the discrete spectrum if a **top” is put on the flow. To examine the source 
of the spurious mode difficulties, suppose that a top is put on the flow at 
some (large) value y = T (say) and the exact outer boundary conditions are 
replaced by the twD conditions (T) = <j>*(T) = 0. In effect, the unbounded 

domain has been replaced by a bounded domain (which could, for exar^le, be 
rescaled to [ 0,l] if desired). It is now well known that stability 
problems on bounded domains have an infinite nunber of discrete modes. 
Consequently, spurious discrete modes are to be expected in global methods 
which impose tops, as the problem they are actually solving does not provide 
a uniform approximation to the true problem for large y, no matter how 
large T is taken. The problem, of course, is that ^ and * in the 
true 0-S problem for boundary layer flows may be exponentially small at T, 
but unless T = (which is not practical computationally), <fr(T) and 

4)'(T) are not identically zero. 

The approach taken in the present work reduces the computational domain 
to 0 < y < Y yet avoids problems with spurious modes by matching the 
computed solution for 0 < y < Y to an appropriate functional form of the 
solution valid asymptotically for y > Y. Specifically, at y * Y, <fr(y) is 
matched to an appropriate exponential solution of equation (3.5) 

E(v.Y) - (3.6) 


wliere, consistent with equation (3.3), for discrete modes 


Y(a,c,R) 


o if Re(a) < Re(X) 
X if Re(X) < Re(o) 


Appropriate choices for y in different situations are: 

V 1) Temporal discrete modes: In this case, a is real, I i 

l.>rge, and a < Re(X). Consequently, the viscous contribution e ^ to 

solutions of (3.3) dies out rapidly relative to the inviscid contribution 
e~^^. For later use in rejecting the possibility of discrete modes close to 
the tanpor^.l continuum, it may be noted that even if 1-cr is small but non- 
zero, |X^ I is still large for large aR, and hence the appropriate choice 
U.7) is always y • a for discrete temporal modes. 



(2) Spacial discrete modes: For these modes, a is complex while the 

frequency w = ac is real. In this case, it is convenient to re-write the 
large parameter as 

ORr:' : 7 rj' 

a 2 . iR(a-o)) + a2 QUALriY ( 3 . 8 ) 

to decide if Y should be a or X . Care is now required in assessing the 
relative speed with which viscous-type and inviscid- type solution of (3.3) 
die out as y -► ® . For fixed values of co and R, Re(a) < Re(X) for most 
values of a . However, because of the complex nature of there are cases 

-A y 

when e dies out slower than e and hence viscous effects penetrate 
relatively far outside the boundary layer. This may be seen by explicitly 
writing a = + i«^ and separating into real and imaginary parts to 

obtain 


X 


1/2 

i(a^R) { 1 - 



i ( 


2ot a . 
r 1 

a.R 

L 



(3.9) 


If R is large while a , and a - m are all order one, then 

o r* X r 

1/2 ■*“^ 1/2 

A » i(a.R) [l - i -Lj ] . (3.10) 

i 

Both real and imaginary parts of [X * i(a - o))/a.]^^^ are order one and 
consequently Re(A ) » 0 ( Oj ^ R )^^2 ^ Hence, e ^ will die out rapidly 
compared to e~°^ and the latter exponential provides the appropriate 
asynptotic form of the eigenfunction for the asymptotic mmerical matching. 
However, as noticed by Houston, Corner, and Ross (1976), if * <•>, to 
lowest order the imaginary part of A is still large, but now 

1/2 

Re(A) 2 o (o./R) « 1 (3.11) 

r i 


As ■ 0(1), the proper asymptotic form of solutions of (3.3) now 
involves e”^^ rather then e"^^. 

A major strength of the present numerical technique is that the outer 
boundary conditions need not be specified apriori and are not intrinsically 
built into the main elements of the code, e.g. into the choice of expansion 
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functions. Rather, boundary conditions enter the present technique through 

separate subroutines. This feature allows specification of the appropriate g 

asymptotic functional form (3.6) of the eigenvalue at Y, but does not ^ 

require that an apriori choice be made for y(a,c,R) {or Rather, 

during the course of the stability calculation, the relative sizes of 

and Re(X ) can be compared, and T can be suitably chosen. As indicated 

above, this flexibility is vital in the case of discrete special modes. 

However, it also allows treatment of continuun modes without the need to 
explicitly change from decaying to merely bounded outer boundaiy conditions. 

This can be done automatically when required based on the computed value of 
Re(X). Mare will be said about continuum modes in the next section. How- 
ever, it is worthwhile emphasizing that the use of asymptotic outer boundary 
conditions not only consistently reduces the computational decision but also 
eliminates the need to treat discrete and continuum modes separately. 

For discrete temporal modes, values of Y from 16 to 18 were found 
adequate to give seven decimal place accuracy in the imaginary part of the 
eigenvalue c. For the special stability problem, corresponding accuracy 
for th«» imaginary part of the eigenvalue a required slightly larger values 
of Y, e.g. 18 to 20. Some typical results of eigenfunctions in both the 
spacial and temporal problems are given in Tables 1 and 2. For the temporal 
problem, the standard test case R « 580 and a » .179 is considered. The 
slightly unstable mode shown has c * .36412 + i 0.0079597. For the spacial 
problem, R » 580 and o) * 0.06520. Table 2 gives the slightly unstable 
mode with a - .17933 - i 0.0033173. 

One further point regarding the present minerical results should be 
considered. Stability calculations are extremently sensitive to errors in 
specifying the mean velocity profile. This is not a problem in the case of 
bounded channel flows or model boundary layer-type flows, such as the asymp** 
totic suction boundary layer profile, which can be specified analytically. 

For the Blasius problem, however, the profile must be obtained nunerically 
as the solution of the boundary value problem 


f’(n) + f(n) f(n) - 0 

f(0) - f'(0) - 0, f* > 1 as n > • 


(3.12) 

ORIGINAL PAbi: 
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Tabl« 1. Un«tabl« dlacrat* taaporal ooda of cha Blaalua atablllty 
problaa for alpha » 0.179 and R ■ 580. Asymptotic boundary 
conditions vara applied at Y • 18, Tbs coaputad value of 
c - 0.36412 + 0.0079597 1. 


X 


FHI 


FHIP 


0*^0000 0 
1*8000 
3»8000 
9*4000 
7*2000 
9*.0000 
10»8000 
12^6000 
14*4000 
16*2000 
18*0000 


*36304E^00 

*4736624>00 

*37839E^00 

*27998E400 

»18069E^00 

*14468E4>00 

*10483i^00 

•79996E-01 

*990342-01 

*39879E-01 


0 * 

•*38399E-02 

*934182-03 

*107992-03 

*208692-09 

*660462—08 

-*366932-10 

*613692-10 

-.292272-11 

-*166412-11 

-•193612-48 


0 * 

*161142«O0 

-*213102-01 

-.636972-01 

-.492802-01 

-.397442-01 

-.298992-01 

-.187692-01 

-*139962-01 

-.989112-02 

-.713772-02 


0 * 

.998212-02 
.199292-03 
-.176712-03 
-* 969922-09 
-.239612-07 
-.369U2-09 
-.394992-10 
.846332-U 
*231312-12 
*139122-39 
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T«bl« 2. Un«t«bl« dlscr«t« spaclal aod* for th« BImIuo ttAblllty 
problem irlth omege ■ 0.06520 end R ■ 58C. Asympeoele boundary 
condlelona var* applied at T ■ 20. The computed value of 
alpha ■ 0.17933 - .0033178 1 eorreapondlng to a wave apeed 
c « 0.36346 * 0.0067246 1. 
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wiiere n is a similarity variable. The form of (3,12) is awkward 
computationally, as it involves an outer condition at infinity. It has been 
known for some time (see, for example, Rosenhead (1963), page 222] that 
(3.12) can be converted to an equivalent initial value problem for a 
function F(s) defined by f(n) * aF(s) where a is a '*homotopy** constant, 
6 = an , and the initial values are F(0) * F; (0) “ 0, F”(0) • 1. Ihe 
homotopy constant a is determined by the limiting behavior of F*(s), i.e. 


a 


^lim F' 


(3.13) 


Bie computation must be carried to sufficiently large s that a is 
determined to the required accuracy. Even slight errors in a may lead to 
a limiting behavior of f'(h) which differs from unity, and this critically 
effects the stability characteristics. 

While relatively simple itself, the system (3.12) provides a cannonical 
example of problems of boundary layer-type which are amenable to homotopy 
methods. It is therefore worth noting that the nunerical determination of 
f(n), f ' (n ) , and f'(n) can be made more efficient through use of a "k- 
homotopy" constant. In particular, define 


f(n ) ■ bF(5 ) with 5 ■ bn 


(3.14) 


Then, F(^ ) satisfies the same equation as f(n), i.e. 

F*»* + FF" ■ 0. 


(3.15) 


Suppose F(C) has initial values 

F(0) - F'(0) - 0 and F''(0) - K. 

Then, (3.14) and the condition F* 1 as n '*' ** give that 

Urn -1/2 

b - ^ F'(e)} and 

f''(0) - K , F'(C)} 
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(3.16) 


(3.17) 



Another quantity of basic interest is the limiting value 0 q 
outflow velocity defined by 


of the Blasius 


I 

f 


, 1 im ^ , 


<5»l5lWAt IT 

<5*^ l»OOR QUALiry 


In terms of the K-variables, this can be expressed as 


For K = 1, the K-variables reduce to the usual horaotopy variables, 
^ 1, tr e relations are 


1/3 

a * K { 


1 im 


c7. 


- 1/2 


(3.18) 


(3.19) 
For K 

(3.20) 


and 


1/3 

8 - K 5 . (3.21) 

To see the advantage in using the R^homotopy technique, let C^(K) 
denote the value of C to which it is necessary to carry the initial value 
calculation in K-variables to achieve sufficient accuracy in the deter- 
mination of the R-constant b. The advantage to choosing R > 1 is that 
is a decreasing function of R. Consequently, by choosing a relatively 
large value of R, the computational domain can be significantly reduced 
for a given accuracy. The relationship betwen R and (R) for eight 

place accuracy is given in Table 3. Corresponding values of f”(0) and 00 
were found to be 

f'(0) - 0.46960020 

a;,d (3.22) 

0Q ♦ 1.2167789 

This value of f"(0), together with f(0) - f'(0) ■ 0, can now be used to 
solve the Blasius equation as a straightforward initial value problem. 
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Table 3. The kappa-homo Copy quantity f (X) for eight place accuracy 
in f"(0). 


X 

f (X) 

1.0 

4.8 

2.0 

3.8 

3.0 

3.5 

5.0 

2.9 

10.0 

2.3 
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THE CONTINUOUS SPECTRUM AND THE POSSIBILITY OF NEARBY DISCRETE 

TEMPORAL MODES 


9 

f 


For conclnuuD modes of Che Orr-Somnerfeld problem, Che ouCer boundary 
ccndicions 0 ac y * are replaced by che weaker condicions 

♦ (y), ♦'(y) bounded as y •. (4.1) 

Grosch and Salwen (1978) have shown Chac on Che conCinuuo Re(X) ■ 0, and 
henro cbe proper funccional form which describes che behavior of ^(y) for 

large y will involve boCh e~^ ^ and e Based on Chis observacion, 

Che Cechnique of replacing Che exacC order condicions (4.1) by asympCocic 
condicions aC Y in nunerical calculaCions based on COLSYS . easily 
exCended Co Che conCinuous specCrun. Indeed, chis can be done in such a way 
chac Chere is no need for apriori specif icacion of whecher a continuous or 
discreCe mode is being compuCed, che decision being me's aucomacically in 
Che course of Che compucacion. Physically, concinutn modes miy be 
identified with free stream vorticity. In particular, the eigenfunctions 
are relatively large outside che boundary layer region. An interesting 
observation, which has not been previously made, is that the appropriate 
parameter X, which is large for the discrete spectrun, is not necessarily 
large on the continuin. For exasple, on the teaporal continuvm, Grosch and 
Salwen have shown that 

c"l-i(l + l^) n/R (4.i.) 

where k is real and non-sero. If this expression is substituted into 

(3.4), the dependence of X on R cancels identically and 

|X| - |k|a (4.3) 

Tints, as k need be merely non-sero, for very lightly daiped continuum 
modes, |X I may in fact be small. 

Miile the temporal eootinuta is simply the straight line ■ 1 and 
< 0, the nature of the special continuum is more complicated due to the 
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complex nature of the eigenvalue o. Corner, Houston, and Ross (1976) have 
found two higher discrete modes of the special problem for the Blasius 
profile which lie fairly close to the special continuum. However, the 
nearby continuous modes are easily distinguished from these discrete modes 
by the behavior of their eigenfunctions. An example of a special continuum 
eigenfunction for the Blasius profile is given in Table 4. For this mode, 

R * 581.1, (I) ■ 0.04649, and the computed eigenvalue is a * 0.047583 + 
0.0085513. The corresponding value of the complex wave speed for this mode 
is c * .94645 -i 0.1701. 

In section 2 of this report, it was proved that there cannot be an 
infinite number of discrete modes in the temporal stability problem for 
boundary layer profiles. Although this is the first theoretical result to 
be proved to date, this fact haa been accepted by most authors for some time 
based on mmerical evidence. Further, most authors feel that for fixed a 
and R, the number of discrete modes is .^all (six or seven for Blasius) . 
Previous nvnerical studies have produced mooos near the continuun, but these 
were usually dismissed as spurious discrete modec attempting to mimic the 
continuun. However, there has been some room for doubt. For example, if 
there were, as perhaps might be inferred from the available nunerical 
evidence, a number of discrete modes lying very close to the continuun it 
appears that neither shooting or expansion techniques could distinguish 
betwen these modes and the continuun functions. A case in point is 
provided by the work of Antar and Benek (1978). They nunerically studied 
the spectrun of Blasius flow by imposing the boundary conditions at some 
finite distance from boundary and attempted to infer the nature of the 
spectrun for the semi~infinite region by systematically increasing the 
distance of the "upper" boundary from the lower boundary. Antar and Benek 
claim to have found two additional groups of discrete modes, which they 
label the P and S families, lying close to the temporal continuun. In 
particular, the P family has close to 1 while the S family was 

initially found to lie along the line C, ■ 0.845. FVirther, the S 
family was conjectured by then to contain an infinite nunber of modes. 

The conjecture of Antar and Benek about the infinite extent of the S 
faiily is disproved by the theoretical results of section 2. the nunerical 
procedure developed in the present work was also used to search for modes 



Table 4. Spaclal continuum mode of the Blaslus stability problem 

for omega » 0.04649 and R * 581.1. Asymptotic boundary conditions 
were applied at Y • 20. The computed eigenvalue Is 
alpha * 0.047583 + 0.0085517 1 corresponding to a wave speed 
c - 0.94645 - 0.17010 1. 
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Table 4. Continued 
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close to the continuum. For a * .179 and R * 580 (in the L scaling) 
values of c attributed to both the P and S fanilies of temporal modes 
by Antar and Benek were used as initial guesses for calculations with 
asymptotic boundary conditions applied at the relatively large value of Y * 
20. The eigenfunctions obtai'ied were all clearly continuim eigenfunctions 
rather than discrete eigenfunctions. A typical calculation is shown in 
Table 5. In light of the present results, it must be concluded that the P 
and S fanilies are spurious discrete modes consequently, there are not large 
nunbers of discrete modes lying close to the continuum. 


BOUNDARY LAYER STABILITY INCLUDING OUTFLOW VELOCITIES 

The parallel flow assumption of linear stability analysis assunes that 
streamwise amplifications of disturbance quantities takes place on a scale 
which is small relative to the scale for boundary la 3 rer growth. Consequent- 
ly, streamwise gradients of mean flow quantities are ignored, e.g. all x 
derivatives of the s t re amf unction ij' for the mean velocity profile. Th.is 
corresponds to neglecting the small outflow (or inflow) velocities of the 
boundary layer and in the case of boundary layer flows over a flat plate 
leads to the usual Orr-Sommerfeld problem. One of the most important flows 
of this type is, of course, the Blasius profile. 

Early theoretical work on the Blasius problem was done by Tollmien 
(1929) and Lim (1945) using piecewise polynomial approximations to the exact 
velocity profile f'(h). When the experimental data of Schubauer and 
Skramstad (1947) first becane available, it was found to be in good agree- 
ment with Che theoretical predictions. Unfortunately, later, more accurate 
nunerical calculations by Jordinson (1971) , who used nunerical solutions 
rather than analytic approximations, for f'(h), showed this apparent 
agreement to be illusory. Attempts were then made to restore agreement by 
taking into account Che non-parallel character of Che basic flow. 

One approach which appeared promising was due Co Barry and Boss (1970). 
They suggest that if i|>(x,y,t) is Che streamfunccion of Che mean flow, Chen 
effects of Che boundary layer outflow velocity can be introduced by retain- 
ing 3t|>^x and neglecting only second and higher streamwise derivatives of 
t. The accuracy of Che mean flow representation in Che stability problem 
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Table 5. Temporal continuum mode of the Blaslus stability problem 
for alpha =* 0.179 and R * 580, Asymptotic boundary conditions 
were applied at Y = 20, Values of c for a supposed discrete 
mode In the P family of Antar and Benek (1978) were used as 
starting values for the calculation. The computed value of 
c - 0.99760 - 0.022995 i. 
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is now the sane as in the Prandtl boundary layer equation. The resulting 

stability problem involves a modified Orr-Soianer feld equation of the form I 

f 

(D^ -o2 ){ (1)2 -a^ ) - RVD}i|; = iaR{ (U-c) (o2-a2) - u}'|» - RV"4» (5.1) 

where the mean profile is now (U,V,0) and V = 0(R~M. 

For the Blasius problem, inclision of the mean boundary layer outflow 
velocity should be destabilizing. However, neutral stability calculations 
of Barry and Ross (1970) based on equation (5.1) led to only a slight reduc- 
tion in the predicted value of R^. Recently, interest in equation (5.1) 
has waned considerably in favor of the multiple-scale approach of Saric and 
Nayfeh (1975). Almost no work has been done on the higher discrete modes of 
the stability problem based on (5.1). 

Recently, Lakin and Reid (1982) have provided a highly accurate 
treatment of the linear boundary layer stability problem based on uniform 
asymptotic expansions involving generalized Airy functions. As a test of 
the as 3 imptotic theory, the stability characteristics of the as 3 mptotic 
suction boundary layer profile, for which 

U(y) “ 1 - e ^ and V*-'/R 0<y<+* (5.2) 

were considered on the basis of both equation (3.1) [v ignored] and equation 
(5.1) [v included]. Results showed that inclusion of the small suction 
component had a significant effect on stability increasing the value of the 

Rc from 47,152 to 54,405. The latter value of R^ was in agreement 
with Hocking's (1975) non-linear treatment of the stability problem 
including suction. 

The approach to stability including V based on equation (5.1) is 
attractive because of its relative simplicity. Based on its competitiveness 
with other approaches in the case of the asymptotic suction boundary layer 
profile, it seemed appropriate to re-examine the stability of the Blasius 
problem based on equation (5.1). This was done using the present numerical 
technique and asymptotic outer boundary conditions involving the appropriate 
large parameter. If outflow is included, the analogue of equation (3.3) as 
y * • is 
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(D?^2)([d . gg /i^Y - = 0 


(5.3) 


where now 

a 2 = iaR(l-c) + o2 + Sq/16 (5.4) 

and Bo is ss (3,18) and (3.22). In particular, the asynrptotic beha^^^*' 
of "viscous- type” solutions of (5.1) as y -► « is now 

g-Xn + ^ (5.5) 

rather than simply e Consequently, the appropriate form of the asynp- 

totic boundary conditions for the discrete modes depends on the relative 
sizes of Be(a) and Re(X-6Q/4) when Re(\) is larger than 

The marginal stability curve computed from equation (5.1) and the dis- 
crete botindary conditions are shown in figure 1. The computed values of the 
critical parameters are Rc » 290.21312 at olq » 0.17 776 and Cq * 

0.40308. The corresponding variation of c with ol along the marginal 
stability curve is given in figure 2. The value of in (5.6) is frac- 

tionally, but not significantly, lower than the value obtained by Barry and 
Ross (1970). For comparison, when the outflow velocity V is neglected, 
Davey (unpublished) has obtained the critical values R * 301.64, a « 
0.17653, and > 0.39664 in the L-scaling. These values are marked with 
an X in figures I and 2. Experimentally observed values of are 261 

[ Schubauer and Skramstad (1947)] and 232 [Boss et al. (1970)]. Figure 1 
only shows the "nose" of the marginal stability curve. However, for values 
of R greater than 500, differences between neutral parameters including 
and neglecting outflow became slight. Further, good agreement is observed 
between the usual linear theory predictions and the experimental dat--^. 

Previous work on the higher modes of the Blasius problem based on equa- 
tion (5.1) [corner, Houston, and Ross (1976)] considered only the special 
case and results were inconclusive. Computations of the higher modes in the 
present work at the standard test values inside the nose of the neutral 
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Variation of c with a along the marginal stability 

curve for the Blaalus problem with proper outflow ln> 

eluded. The circle denotes the values at R . The x 
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marks the corresponding point when outflow is neglected 



curve did not suffer from convergence problems and showed very little dif- 
ference in decay rates when V was retained or neglected. The eigenfunc- 
tions of mode 3 of the temporal stability problem for a = 0.179 and R = 

580 are given as an example in Table 6. With V ignored the computed 
eigenvalue is c * 0.48393 - i. 19207 while with V included the eigenvalue 
is c = 0.48698 - i. 19208. 

The present results show that while linear stability theory is excel- 
lent for predicting general stability characteristics, the inclusion of the 
velocity V along the lines of equation (5.1) is not sufficient to recon- 
cile theory and experimental results for boundary layer profiles involving 
outflow. When contrasted with the success of equation (5.1) in the case of 
the asymptotic suction boundary layer profile, a tentative conclusion may be 
drawn: For boundary layer flows where V is destabilizing, boundary layer 

growth effects may require the use of more complicated multiple-scale 
approaches to refine linear stability predictions. However, if V is stab- 
ilizing, good results can be expected from the linear stability problem 
based on equation (3*1). In particular, if suction is used on a flat plate 
(back from the leading edge) to retard instabilities and retain a laminar 
flow, the resulting local V velocity will be stabilizing and equation 
(5.1) should give accurate stability predictions. 
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Table 6. The effect of partially relaxing the parallel flow 
assumption on mode 3 in the discrete temporal spectrum 
of the Blasius stability problem with alpha * 0.179 and 
R “ 580. Asymptotic boundary condition was applied at 
Y - 18. 


A. The eigenfunction based on equation (3.1) where V is 

ignored. The computed eigenvalue is c - 0.48393 - 0.19207 i. 
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Table 6. Continued. 
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3. The eigenfunction based on equation (5.1) where V Is Included. 
The computed eigenvalue is c ■ 0.48698 - 0.19208 1. 
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